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The discovery of high-energy (TeV-PeV) neutrinos from gamma-ray bursts (GRBs) would shed light on their 
intrinsic microphysics by confirming hadronic acceleration in the relativistic jet; possibly revealing an accel- 
eration mechanism for the highest energy cosmic rays. We describe an analysis featuring three models based 
upon confronting the fireball phenomenology with ground-based and satellite observations of GRB030329, 
which triggered the High Energy Transient Explorer (HETE-II). Contrary to previous diffuse searches, the 
expected discrete muon neutrino energy spectra for models 1 and 2, based upon an isotropic and beamed emis- 
sion geometry, respectively, are directly derived from the fireball description of the prompt 7-ray photon energy 
spectrum, whose spectral fit parameters are characterized by the Band function, and the spectroscopically ob- 
served redshift, based upon the associated optical transient (OT) afterglow. For comparison, we also consider a 
model (3) based upon averaged burst parameters and isotropic emission. Strict spatial and temporal constraints 
(based upon electromagnetic observations), in conjunction with a single, robust selection criterion (optimized 
for discovery) have been leveraged to realize a nearly background-free search, with nominal signal loss, using 
archived data from the Antarctic Muon and Neutrino Detector Array (AMANDA-II). Our preliminary results 
are consistent with a null signal detection, with a peak muon neutrino effective area of 80 m^ at ^ 2 Pe V and 
a flux upper limit of 0.150 GeV/cm^/s for model 1. Predictions for IceCube, AMANDA'S kilometer scale 
successor, are compared with those found in the literature. Implications for correlative searches are discussed. 

1. Introduction 

Neutrino astronomy may provide us with a new glimpse at the internal processes of gamma-ray bursts (GRBs). 
The Antarctic Muon and Neutrino Detector Array (AMANDA), which has been calibrated upon atmospheric 
neutrinos, has demonstrated the viability of high energy neutrino astronomy by using the ice at the geographic 
South Pole as a Cherenkov medium. Canonical fireball phenomenology 1 1 1, in the context of hadronic accel- 
eration, predicts correlated MeV-EeV neutrinos from GRBs via various hadronic interactions j2|. Ideal for 
detection are TeV-PeV neutrinos in coincidence with prompt 7-ray emission, resulting in a nearly background- 
free search. A detailed description of the modeling techniques and an ongoing analysis featuring correlated 
neutrino searches of individual GRBs from the Burst and Transient Source Experiment (BATSE) are de- 
scribed elsewhere f3', "^l. Here, we report on a complementary search for correlated leptonic (j/^, i?^) emis- 
sion, using models based upon the unique discrete electromagnetic characteristics and emission geometry of 
GRB030329, gleaned directly from satellite and ground-based observations. This represents a novel departure 
from searches |4l|5]|6| based upon a diffuse formulation |7 1, which utilize averaged burst parameters. 

2. GRB030329: Electromagnetic Emission & The GRB-Neutrino Connection 

On March 29, 2003, at 11^37" 14.^67 (UTC), HETE-II was triggered by GRB030329 (H2652), a watershed 
transient which confirmed the connection between a core collapse Type Ic supernova and long duration GRB. 

"Adapted from a contribution to the Proceedings of the 29''' International Cosmic Ray Conference in Pune, India (August 2005). 
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Electromagnetic investigations of the prompt 7-ray and multi-wavelength afterglow emissions associated with 
GRB030329 abound in the literature (see tabled, making it a perfect specimen for study. Via forward folding 
deconvolution, the photon energy spectrum was fit to the Band function (Sj, an empirically derived power law 
with smooth transition. For spectral indices a > —2 and (3 < —2, the characteristic peak energy is defined as 
Ep = [(2 + Q!)e!j,] (a — /3)~^, where is known as the photon break energy 1 8 1. Hence, using standard error 
propagation, we find that e^J, = 115.6 ± 9.9 keV. Doppler spectroscopy of the OT afterglow revealed a redshift 
measurement, which, under an assumed Acdm cosmology,' placed GRB030329 at a luminosity distance of 
2.44j]Q^g X 10^^ cm. Coupled with the peak energy flux, this implies an intrinsic peak isotropic luminosity of 
j^iso ^ 5.241q'77 X 10^° ergs /s in the 30-400 keV energy band pass. Evidence for anisotropic emission, in the 
form of a two component break in the afterglow spectrum, was revealed by radio calorimetry and is consistent 
with collimated prompt emission within a jet of opening half angle 6 jet- This requires a beaming fraction 
correction, which reduces the intrinsic peak luminosity to L^'^* = Ll^° {1— cos 9 jet ) ~ 1.99^Q'2g x IC^^ergs/s. 
Extended calorimetry provided an estimate for the fractions of shock energy imparted to the electrons (eg) and 
magnetic field (es). Table^summarizes the observed electromagnetic properties used in this analysis. 



Table 1. Electromagnetic Properties of GRB030329: Prompt 7-ray and Multi-wavelength Afterglow Emission 



Paraineter(s) 


Value(s) 


Reference 


Positional Localization (aj2ooo, <5j2ooo, c^h) 
Trigger (T) & Duration (T90) [30-400 keV] 

Energy Fluence (F^) [2-400 keVl 
Band Parameters {a, 13, Ep) [2-400 keV] 
Peak Energy Flux (*!^°'*) [30-400 keVl 
Spectroscopic Redshift (z) 
Jet Opening Half Angle (6 jet ) 
Election & Magnetic Field Energy Fractions 


161.2081646°, 21.5215106°, 3.0 x 10"^ ° 
41, 834.67 SOD, 22.8 ± 0.5 s 
l.Q30t°o°o\t X 10-* ergs/cm^ 
-1.32 ± 0.02,-2.44 ± 0.08, 70.2 ± 2.3 
~ 7 X 10"'' ergs/cm^/s 
0.168541 ± 0.000004 

~ 5° f» 0.09 rad 
Ce ~ 0.19, CB =s 0.042 


Taylor et al., GCN Report 2129 
Vanderspeket al., ApJ 617: 125I-I257 (2004) 

Sakamoto et al., astro-ph/0409 1 28 
Vanderspeket al., ApJ 617: 1251-1257 (2004) 
Vanderspeket al., GCN Report 1997 
Bloom et al., GCN Report 2212 
Berger et al.. Nature 426, 154-157 (2003) 
Frail et al., ApJ 619, 994-998 (2005) 



The generic mechanism responsible for the super-Eddington luminosities associated with GRBs is the dissipa- 
tion, via shocks, of highly relativistic kinetic energy, acquired by electrons and positrons Fermi accelerated in 
an optically thick, relativistically expanding plasma, commonly referred to as a. fireball. The acceleration of 
electrons in the intense magnetic field of the fireball leads to the emission of prompt non-thermal 7-rays via 
synchrotron radiation. The temporal variability (t^, ^ 10 ms) in the fight curves imply compact sources. In 
order to ensure a transparent optical depth to photons of energy g™"^ k, ICQ MeV, a minimum bulk Lorentz 
boost factor (F) was assigned (see equation l2dt. Hadronic acceleration within the ambient photon field pro- 
duces TeV-PeV leptons, via the following photomeson interaction: 



p + 7 A+ ^ 7r+ + [n] — > ^+ + v^, + i^e + Vf_L + v^^ (1) 

Hence, these neutrinos are expected to be spatially and temporally correlated with prompt 7-ray emission. The 



^AcDM cosmology: Ho = 72 ± 5 kiWMpc/s, n,„ = 0.29 ± 0.07, = 0.73 ± 0.07 [Spergel et al., ApJS 148, 175-194 (2003)]. 
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parameterization of the neutrino energy spectra {dNi, / de^, = <i>,^ ) traces the prompt photon energy spectra, 
as illustrated in figure[l] and is defined as follows ^3,,9J: 



X < 







-I3~l 

-Q-1 



4 < < 4 



(2a) 



~ w/n^^ ~ 9.86 X lO"'^ GeV/cmVs (2b) 



8e, ln(10)Tc 



90 



ri5^i',-24.Mev(l + ^) 



A - 0-2 X ^ J'"' , ~ 0.77 (2c) 



r > 276 [L^,52i;i2CMey(l + z)J ' « 178 (2d) 



7x10^ 



(1 + ^) 



2 ,6 



'^7,MeV 



GeV w 1.404951 x 10'' GeV 



(2e) 



GeV « 7.9832941 x 10^ GeV 



(2f) 



The values in equations I2bl2t1 are given for model 1, where = L^,52 • 10^^ ergs/s, T = T2.5 ■ 10 



2.5 



tu = tu,_2 • lOms, e!^ = eJ^ Mev • IMeV, and 
7-ray photon observables. 



100 MeV. Note the explicit dependence on discrete 



3. Neutrino Astronomy with AMANDA-II: Analysis, Results & Discussion 

AMANDA-II is comprised of 677 optical modules buried at depths between 1500-2000 m. The background 
consists of cosmic ray induced down-going atmospheric muons, detected at a rate of ^ 100 Hz, with a pertur- 
bation of atmospheric neutrinos, detected at a rate of ^ 10~* Hz. The astrophysical neutrino signal, detected 
via charged current interactions such as: j/^ + iV ^ fi^ + X,is isolated by utilizing topologically identified 
up-going muon events, which are reconstructed by a maximum likelihood method. On-source, off-time data 
were used to estimate the stability of the background rate in order to maintain blindness, facilitating an un- 
biased analysis. After data filtering (see | lOj), the total off-time background, excluding a 10 minute blinded 
window centered upon the trigger time, was consistent with a Gaussian fit and accrued 24, 972 ± 158 events 
over a 57,328.04 second interval, resulting in a rate of 0.436 ± 0.003 Hz. Based upon a visual inspection of the 
light curve, a conservative search window of 40 seconds (beginning at T) was chosen. Hence, 17.44 ± 0.12 
background events (rifc) were expected on-time in AMANDA-II prior to any quality selection. Signal neutrino 
spectra were simulated for three models (see figure^ by propagating a total of 440, 000 muon neutrinos 
(i/^ + Df^) from an error box in the sky defined by the spatial localization of the radio afterglow (see tableQ. 
Event quality selection was optimized for the best limit setting and discovery potential by minimizing the 
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model rejection factor (MRF) 1 1 1 1 and the model discovery factor (MDF) 1 12 1, respectively. Although multi- 
ple observables were investigated, a single criterion emerged, based upon the maximum size of the search bin 
radius (^f), i.e. the space angle between the reconstructed muon trajectory and the GRB's position. 

Our search, optimized for 5cr discovery (requiring 4 events within 5* < 11.3°), is consistent with a null signal. 
Upper limits, summarized with our results in table |2] do not constrain the models tested using AMANDA-II. 
Effective neutrino and muon areas are given in figure^ The number of expected signal events in IceCube {Ng) 
for model 1 is consistent with 1 13 1, when neutrino oscillations are considered. For model 3, our results for Ng 
are in agreement with |9 1 and 1 14|, when one adjusts for the assumptions of |7 1. Selection based upon ^ was 
robust across the models tested in AMANDA-II, as illustrated in figure|2] However, the MRF/MDF and hence 
the limit setting/discovery potential was strongly model dependent, varying by over an order of magnitude for 
models 1 and 3. Furthermore, using the same theoretical framework, the response of AMANDA-II and IceCube 
to spectra based upon discrete and average parameters are discrepant in mean neutrino energy and event rate by 
over an order of magnitude as illustrated in figure and table |2l Such variance in detector response unequiv- 
ocally demonstrates the value of a discrete modeling approach when making correlative neutrino observations 
of individual GRBs, especially in the context of inferred astrophysical constraints, in agreement with |9J. 



Table 2. Summary of Results for GRB030329. Primes indicate value after selection for IceCube (N^,) and AMANDA-II 
{vix). Superscripts indicate A=MRF and B=MDF optimization. Average upper limit (sensitivity) and flux upper limit values 
are given for AMANDA-II. The effects of neutrino flavor oscillations have been included. 



Neutrino 


Search 


Background 


Expected Number 


Observed 


Optimization 


Average 


Flux 


Flux 


Bin 


Expected in 


of Neutrino 


Events in 


Method for 


Upper limit^ 


Upper 


Model 


Radius 


AlVlANDA-lI 


(Signal) Events 


AMANDA-11 


AMANDA-11 


[Sensitivity] 


Limit 






rib, n^' ,n^' 


Ns, n^, nf 


nabs , nfb's 


MRF, MDF 


GeV 


Gey 


1 


21.3, 11.3 


17.44, 0.23, 0.06 


0.1308,0.0202, 0.0156 


15,0 


152,424 


0.157 


0.150 


2 


18.8, 11.3 


17.44, 0.17, 0.06 


0.0691,0.0116, 0.0092 


15,0 


256,716 


0.041 


0.039 


3 


18.5, 11.3 


17.44, 0.17, 0.06 


0.0038, 0.0008, 0.0006 


15,0 


3864, 10794 


0.036 


0.035 
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Figure 1. Upper panel - Prompt neutrino energy flux for models 1 (solid black), 2 (dashed blue) and 3 (dot-dashed red), 
based upon equation |23 for GRB030329. Lower panel - Detector response for models 1 (black), 2 (blue) and 3 (red) 
for AMANDA-II (solid) and IceCube (dashed). Note how the variance in the neutrino spectra (upper panel), due to the 
variance of electromagnetic parameters used in equation|2^ directly translate into the variance in detector response (lower 
panel), which exceed an order of magnitude, as manifested in observables such as the expected mean neutrino energy and 
rate. The effects of neutrino flavor oscillations have been included. 
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Figure 2. AMANDA-II signal efficiency/background rejection for GRB030329 models 1-3 using MRF (vertical dashed 
blue, dashed-dot green and dashed-dot-dot red) and MDF (vertical black dotted) optimizations. Both MRF and MDF 
selections reject ~ 99% of the background while retaining ~ 86% and ~ 77% of the signal, respectively. 
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Figure 3. Upper panel - The effective muon neutrino area. A peak of ~ 80 and ~ 700 occurs at ~ 2 PeV 
for AMANDA-II and IceCube, respectively. Lower panel - The effective muon area for energy at closest approach to the 
detector. At ^ 200 TeV, the effective muon area for IceCube, ~ 1 km^, exceeds that of AMANDA-II, ~ 100, 000 m^ by 
an order of magnitude. For both panels, MDF optimized AMANDA-II results for models 1 (dashed black), 2 (dashed blue) 
and 3 (dot-dashed red) and predicted IceCube (solid black) curves are illustrated for GRB030329 with fcooo ~ 22°. 



